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We present a perturbative method for calculating the optical Hall conductivity in a tight-binding 
framework based on the Kubo formalism. The method involves diagonalization only of the Hamil- 
tonian in absence of the magnetic field, and thus avoids the computational problems usually arising 
due to the huge magnetic unit cells required to maintain translational invariance in presence of a 
Peierls phase. A recipe for applying the method to numerical calculations of the magneto-optical 
response is presented. We apply the formalism to the case of ordinary and gapped graphene in a 
next-nearest neighbour tight-binding model as well as graphene antidot lattices. In both case, we 
find unique signatures in the Hall response, that are not captured in continuum (Dirac) approxima- 
tions. These include a non-zero optical Hall conductivity even when the chemical potential is at the 
Dirac point energy. Numerical results suggest that this effect should be measurable in experiments. 

PACS numbers: 81.05.ue, 78.20.Ls,78.20.Bh 



I. INTRODUCTION 

Since the experimental discovery of graphene^i the 
honeycomb lattice has been the subject of intense 
researchi^^— Graphene displays unique properties in an 
external magnetic field, with a non-equidistant Landau 
level structure and a zeroth Landau level energy which is 
independent of the magnetic field strength.—"— The Lan- 
dau level structure is reflected in the half-integer quan- 
tum Hall effect observed in graphene?^ which, due to 
graphene's large cyclotron gap, has been observed at 
room temperature 1^ These and other remarkable features 
of graphene emerge quite naturally from a low-energy, 
continuum description of graphene, the so-called Dirac 
approximation, which is based on a linearization of a 
nearest-neighbor (NN) tight-binding (TB) model near 
the high-symmetry _ftr-pointsJ^ 

While many properties of graphene are correctly de- 
scribed by the Dirac model, it nevertheless fails in cer- 
tain respects. A simple example is found in the optical 
response displaying a clear saddle point resonance around 
4.4 eV that is absent in a linearized mode h^^d^ More sub- 
tle effects such as an orbital magnetic susceptibility away 
from the Dirac point have been identified as consequences 
of lattice effects lost in a continuum approach Here, we 
demonstrate that under certain circumstances the Hall 
conductivity, which is routinely applied as an important 
tool to characterize graphene experimentally ;^' ^^d^ can- 
not be accurately described by the Dirac model. Specif- 
ically, the perfect electron-hole symmetry of the Dirac 
model, present also in the non-linearized NN TB model, 
results in an optical Hall conductivity which is identi- 
cally zero unless electron-hole symmetry is broken by 
moving the chemical potential away from the Dirac point 
energyji^iii In the present work, we predict that going be- 
yond the Dirac model by including next-nearest neighbor 



coupling in a full TB model yields an appreciable optical 
Hall response even when the chemical potential coincides 
with the Dirac point. This result is demonstrated using 
a novel perturbative technique that allows us to evaluate 
the magneto-optical response in an atomistic model for 
arbitrarily small field strengths. Thus, we identify a sig- 
nificant lattice effect, for which the Dirac model predicts 
a null result. We note that, interestingly, this deviation 
occurs for energies well within the range of the linearized 
band structure of graphene. It is thus not a result of sim- 
ply probing the band structure beyond the validity of the 
linearized model, as is the saddle point resonance men- 
tioned above, but is rather a strong signature of broken 
electron-hole symmetry. 

To go beyond the Dirac model, we return to the 
tight-binding model from which the Dirac approxima- 
tion emerges. The effect of a magnetic field can then 
be included via a Peierls substitution. The trouble with 
this method is that the periodicity of the Peierls phase, 
for realistic magnetic field strengths, is usually orders of 
magnitude larger than the lattice constant. Thus, cal- 
culations must be made on magnetic unit cells hundreds 
or thousands of times larger than the Wigner-Seitz cell. 
In the case of bulk materials, where the unit cell con- 
sists of just a few atoms, this problem may be over- 
come for large, but reasonably realistic magnitudes of the 
magnetic field. However, for nanostructured graphene 
materials, where the Wigner-Seitz cell may itself con- 
tain hundreds of atoms, direct diagonalization of the re- 
sulting Hamiltonian is not feasible. Several numerical 
methods have previously been suggested to overcome this 
problem.—"— However, all these methods eventually fail 
at sufficiently small magnetic fields, because of the diver- 
gence of the size of the magnetic unit cell. 

In this paper, we present a perturbative approach to 
calculating the optical Hall conductivity in TB models 
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on a honeycomb lattice. The approach requires diag- 
onahzation only of the Hamiltonian in absence of the 
magnetic field, and thus circumvents the problem associ- 
ated with the periodicity of the Peierls phase. We apply 
the formalism first to ordinary and gapped graphene, il- 
lustrating clear and qualitative deviations from a Dirac 
approximation, as discussed above. Finally, to illustrate 
the power of the perturbativc formalism, we apply the 
method to an example of nanostructured graphene, in 
this case graphene antidot lattices j ^-'^i^^ and once again 
find qualitative differences in the optical Hall conductiv- 
ity compared to a simple continuum treatment. 



II. PERTURBATIVE METHOD 
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FIG. 1: (Color online) Unit cells used for the analytical and 
numerical calculations of the optical Hall conductivity of (a) 
bluk graphene, and (b) a {4, 2} graphene antidot lattice. The 
grey shading indicates the size of the unit cell, with high- 
lighted carbon atoms included in the unit cell. Note that the 
rectangular unit cell for the GAL is chosen for computational 
convenience only. 



The derivation of the main result of the perturbative 
treatment is given in full details in Ref. [231 The ap- 
proach is based on the strategy developed in Refs. \24 
and [2^ The starting point is a TB approximation of the 
honeycomb lattice without magnetic field. The optical 
conductivity is then evaluated using the Kubo formal- 
ism, with the effect of the magnetic field included per- 
turbatively via a Peierls substitution. That is, a phase 
is added to the hopping terms Uj between atomic sites 
i and j, such that tij — > tije'^'^, with the phase given as 

4> = e/h X J^;' A • dl. Here, Ri and Rj denote the posi- 
tions of the atomic sites, while A is the magnetic vector 
potential. Including the effect of the magnetic field to 
first order, the result for the optical Hall conductivity 
reads as 



X Ti [T,y{z) + T,y{z + fi^)]}. 



(1) 



Also, we note that all even powers of the expansion in the 
magnetic field strength are zero, so the equation is valid 
up to third order in the magnetic field strength. Here, B 
is the magnetic field strength, e the electron charge, Hll) 
is the photon energy and f{z) is the Fermi-Dirac distri- 
bution function. The first integral is over the Brillouin 
zone, while the last integral is to be taken along a con- 
tour C, which should enclose the entire energy spectrum 
of the Hamiltonian, H, while C ± fku should not contain 
the spectrum. Here, and in what follows, we include a 
small imaginary part in huj = huo + to account for 
broadening. Note that a spin degeneracy factor has not 
been included, and the final trace thus includes tracing 
over spin degrees of freedom. This trace is over the op- 



erators T^y{z) = T^l\z) + T^t)'{z), with 
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Here, we have defined the derivatives of the Hamiltonian 
Hi = dH/dki and Hij = d^H/dkidkj, with similar def- 
initions for the Green's functions G ~ {H ~ z)~^ and 
G = {H — z + fujj)^^ . We stress that because the mag- 
netic field is treated as a perturbation, the Hamiltonian 
appearing in these expressions is the Hamiltonian in ab- 
sence of the magnetic field. For numerical calculations, 
this is a significant advantage of this method, as we will 
demonstrate in more detail below. 



III. LINEARIZED GRAPHENE 

Before turning to the full TB model, we wish to apply 
the above approach to demonstrate that, indeed, the Hall 
conductivity in the linearized model vanishes if electron- 
hole symmetry remains unbroken. We consider the TB 
Hamiltonian corresponding to the unit cell of graphene 
shown in Fig. [IJi. For generality, in addition to ordi- 
nary bulk graphene, we consider also a gapped graphene 
model, where a staggered mass term ±A is added to the 
on-site energies, with the sign alternating between the 
two sublattices of graphene. Pristine graphene may ex- 
hibit a gap due to excitonic effects i^^i^^ However, here we 
will focus on magnitudes of the gap that are relevant for 
nanoengineered graphene, such as, e.g., graphene antidot 
latticeSi ^^'^^ We stress that the results obtained remain 
qualitatively the same for any magnitude of the gap. Lin- 
earized around the /^-point, we find that the optical Hall 
conductivity can be evaluated as 
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where we have introduced the zero-frequency graphene 
conductivity CTq = e^/Ah and the cyclotron frequency 
Wc = V F \J "^^eB / with vp the Fermi velocity. We note 
that the K and K' valleys contribute equally, and that 
a factor of 2 to account for this valley degeneracy has 
already been included in the above equation. This result 
agrees with the low field strength limit of previous analyt- 
ical results derived by Gusynin et al.i^ The full details 
of the derivation of Eq. Q are given in the appendix. 
We note that the final term /'(— e) — /'(e) oc sinh(/i/fcr) 
with ^ the chemical potential while kT is the thermal 
energy. This demonstrates how the optical Hall conduc- 
tivity is identically zero in the symmetrical case, where 
the chemical potential sits at the Dirac point energy. Be- 
low, we will demonstrate how this zero-result is drasti- 
cally altered when going beyond the continuum (Dirac) 
treatment of graphene. 



IV. NUMERICAL RESULTS 

The analytical result for linearized graphene, presented 
above, is interesting in its own right, and serves as a way 
of validating the perturbative approach. However, the 
real power of the method lies in the fact that because the 
expression in Eq. ([l} is given in terms of the Hamiltonian 
without magnetic field, numerical TB calculations can be 
performed on a drastically smaller unit cell than using 
the standard method of Peierls substitution in a non- 
perturbative manner. Peierls substitution necessitates a 
magnetic unit cell large enough to ensure periodicity of 
the magnetic phase factor added to the hopping terms. 



For graphene, this leads to a scaling of the total number 
of carbon atoms as TV ~ 316 • 10'^ T x , rendering 
realistic magnetic fields quite difficult to manage using 
this methodjii 



A. Numerical recipe 

To arrive at an expression suitable for numerical sim- 
ulations, we first note the identity Gi = —GHiG, with 
Hi = dH/dki, and a similar definition for the Green's 
functions G = [H — z)"^ and G = {H — z + hoj)"^ . Using 
this identity we write the trace in the eigenstate basis as 
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where we have introduced z± — z ±fvjj and 



Ml^pg - {mmn) {n\H,\p) {p\Hk\q) {q\Hi\m) . (6) 
We can now perform the contour integration using the 
residue theorem. To ease notation we define Emn = Em — 
i5„, Af^i^^f^"''^''' = M^a^'pq-M^;^'/ aswellas J„„ = 1- 
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delta. We then arrive at the rather lengthy expression 
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where we have introduced fl = huj. In a similar fashion, 
we write the trace over the second operator as 
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The residue theorem then leads to 

dz = 
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The trace of the operators with the shifted argument is 
obtained in a similar fashion. By substituting z ^ 
and z_ — > z in Eq. ([5|) and relabeling slightly, one can 
show that the contour integral f{z)Tr[Txl\z + fl)] is 
given by Eq. ([7]) , if one substitutes 



l\/fxxyy,yxxy 
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mnpq 
ji^xyxy,xxyy 

mnpq 

.(2) 
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Similarly, §^ f[z)Tv[Tiy\z + n)] is given Eq. (dU]), with 
the substitutions 
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(12) 



In this way we have arrived at expressions for the con- 
tour integral of the trace in Eq. (1) in terms of sums 
over the eigenstates of the Hamiltonian without mag- 
netic field. These sums can be truncated to include only 
states near the Fermi energy. For small magnetic fields, 
this method is drastically faster than direct diagonaliza- 
tion of the Hamiltonian with magnetic field, the size of 
which diverges as the magnetic field strength is reduced. 
In all numerical results presented below, we set the ther- 
mal energy to kT = 0.025 eV and include a broadening 
of hV = 0.05 eV. 



B. Graphene 

We now consider a next-nearest neighbor (NNN) TB 
model of gapped graphene, defined via the Hamiltonian 



iJ(k) 



t'g(k) + A -i/(k) 
-tr(k) ^'.9(k)-A 



(13) 



parametrized by the nearest and next-nearest neigh- 
bor hopping parameters t = 3 eV and t' = 
0.3 eV, respectively. Here, /(k) = g*fea=ac _|_ 
2e-'''-°-/2cos(\/3fcyac/2), while ^(k) = 2 cos(%/3fcyac) + 
4 cos(3fca;ac/2) cos(\/3fcj,ac/2), where = a/-\/3 is the 
carbon-carbon distance. While the hopping parameters 
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FIG. 2: (Color online) Optical Hall conductivity, (Jxy, as a 
function of photon energy for ordinary graphene and gapped 
graphene with increasing values of the mass term. The chem- 
ical potential is in the middle of the gap. The conductivity is 
shown relative to the zero- frequency conductivity of graphene, 
(70, times the square of the cyclotron energy, fkiJc- The solid 
(dashed) lines show the real (imaginary) part of the conduc- 
tivity. To ease visibility, the response of ordinary graphene 
has been cropped. The DC value of ordinary graphene is 
cTxyi'd) — {hjJc)^ao X 1.03 eV~^. The thin line shows the real 
part of the Hall conductivity for A = 0.1 eV and a magnetic 
field of B = 26.3 T. 



can vary slightly depending on which ab initio results 
they are fitted to, we note that the exact value of the 
hopping terms do not alter out results qualitatively. We 
set the on-site energy to zero. For t' = 0, this model has 
electron-hole symmetry, which is inherited in the Dirac 
model discussed above. We note that linearization of any 
TB model will inevitably result in perfect electron-hole 
symmetry. As discussed above, in the fully symmetrical 
situation, where the chemical potential sits at the Dirac 
point energy, the off-diagonal conductivity is identically 
zero for any such model. This can be proven on quite gen- 
eral terms for all TB models exhibiting tt-tt* symmetry, 
for which contributions from conjugate transitions ex- 
actly cancel in the fully symmetrical caseJ^iii Introduc- 
ing next-nearest neighbor coupling breaks electron-hole 
symmetry and, as we will now demonstrate, leads to a 
markedly different magneto-optical response of graphene. 

In Fig. [5] we show the calculated optical Hall con- 
ductivity of graphene, with the chemical potential at 
the Dirac point energy. While the Dirac approxima- 
tion (and nearest-neighbor TB) predicts a zero response 
in this case, our NNN TB model suggests a clear reso- 
nance at hu = 2A. This drastic deviation from the Dirac 
model is due to the broken electron-hole symmetry, which 
means that conjugate transitions in the optical response 
no longer cancel entirely The strength of the resonance 
decreases as the magnitude of the band gap is increased, 
in agreement with previous results showing that a suf- 
ficiently large band gap effectively quenches the effect 
of the magnetic field, provided A ^ fkUc'^ However, 
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FIG. 3: (Color online) Optical Hall conductivity, a^y, as a 
function of photon energy for ordinary graphene and gapped 
graphene with increasing values of the mass term. The chem- 
ical potential is in the middle of the gap. The conductivity is 
shown relative to the zero- frequency conductivity of graphene, 
fjo, times the square of the cyclotron energy, fiajc. The solid 
(dashed) lines show the real (imaginary) part of the conduc- 
tivity. The thin lines show the corresponding results obtained 
via the semi-analytical expression derived in the main text. 



the magnitude of this correction to the Dirac response is 
quite appreciable, suggesting that these deviations from 
the Dirac model should be measurable in experiments. 
We note that, as expected, numerical calculations show 
similar results for a nearest-neighbor model if overlap be- 
tween neighboring tt orbitals is included. 

For comparison with the perturbativc results, we also 
show the Hall conductivity for A = 0.1 eV and a mag- 
netic field strength oi B = 26.3 T, calculated using 



standard, non-perturbative tight-binding methodsJi We 
note that these calculations involve a Hamiltonian with 
12000 X 12000 elements for what is a quite strong mag- 
netic field. The relationship with the perturbativc re- 
sult is evident in the figure, and illustrates the fact that 
the perturbativc results still have predictive power for 
the strength of the response even in substantial magnetic 
fields. In particular, the perturbativc results correspond 
to an averaging of the oscillations occurring due to in- 
dividual Landau levels, which for smaller magnetic field 
strengths could presumably be caused by a broadening 
of the order of the cyclotron energy. 

To further corroborate these findings, we will now de- 
rive an approximate, semi-analytical expression for the 
optical Hall conductivity in the next-nearest neighbor 
model. We first note that linearizing the NNN Hamil- 
tonian in Eq. ((T5)) near the K point yields the same re- 
sult as the nearest-neighbor model, except for a constant 
diagonal term. Instead, we proceed by expanding the 
diagonal NNN term to second order near the K point, 
yielding the approximate Hamiltonian 
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iKy) 

-A 
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where we have defined Ki = tkiac and introduced the 
2t' /t^, quantifying the perturbation due 



The eigenvalues of this Hamiltonian 
We now proceed in a 
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manner similar to that of the appendix, i.e., we evaluate 
the trace of Txy{z), integrate out the angular component 
of the Brillouin zone integral and then use the residue 
theorem to perform the contour integral over z in Eq. ([Ij . 
In this manner, we find that the optical Hall conductivity 
in the NNN model is approximately given by 



Oxy\^) = Co 5- 



2 fco 1, 



r(3e2 _ A2)n2^2 [^(^^) _ ^(^_)] + a^{e)r{E^) + a+(e)/'(i?+) 



A 



e2(;i2^2 _4£2) 

where we have introduced 

a±(e) = e {2e(Te ± 1) (l - (re)^) (e^ - A^) - h^u? (±e - \t [(re)^ - 1 ± 2Te 

I 



£2 _ a: 



(15) 



(16) 



where E± ^ ^T{e^ ~ A^) ± e, with e = ^ + |k2. 

In Fig. [3] we compare the numerical results with those 
obtained by numerical integration of the analytical result 
derived above. We note that there is excellent agreement 
between the two methods. 



C. Graphene antidot lattices 

Nanostructured graphene systems, with Wigner-Seitz 
unit cells containing on the order of hundreds of atoms, 
are practically impossible to handle using direct diago- 
nalization of the TB Hamiltonian in the presence of a 
realistic magnetic field. To illustrate the power of the 
perturbativc method presented above, we now consider 
the magneto-optical response of periodically perforated 
graphene, so-called graphene antidot lattices (GALs)<^ 
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FIG. 4: (Color online) Optical Hall conductivity, axy, as a 
function of photon energy for the {4, 2} graphene antidot lat- 
tice. The conductivity is shown relative to the zero- frequency 
conductivity of graphene, ctq, times the square of the cyclotron 
energy, fkuc- The blue (red) lines show the real (imaginary) 
part of the conductivity. The thick lines show the results of 
the perturbative method applied to the GAL structure, while 
the thin lines are results of a gapped graphene Dirac model 
with a band gap corresponding to the GAL. 



The low-energy spectrum of these structures can be quite 
accurately described in a gapped graphene model, by fit- 
ting the mass term to coincide with half the magnitude 
of the GAL band gapi^ We now compare TB results to 
such a continuum description of GALs. For these results 
we ignore next-nearest neighbor coupling, to illustrate 
how deviations from a continuum approximation emerge 
even in the simplest nearest-neighbor TB treatment. We 
denote a given GAL structure as {L,R}, where L is the 
side length of the hexagonal Wigner-Seitz cell, while R 
denotes the radius of the circular hole in the middle of 
the cell, both in units of the graphene lattice constant. 
We consider a geometry for which the superlattice basis 
vectors are parallel to the carbon-carbon bonds, as such 
structures always exhibit band gaps<^ As an example. 
Fig. [TJd shows the computational cell of a {4, 2} GAL, 
highlighted with gray shading. We note that the rect- 
angular unit cell contains 144 carbon atoms. For com- 
parison, a standard non-perturbative calculation of the 
magneto-optical properties would require a magnetic unit 
cell consisting of 72000 carbon atoms, even for a substan- 
tial magnetic field strength of 40 T. 

In Fig. U we show the optical Hall conductivity of 
the {4, 2} GAL calculated using the perturbative ap- 
proach. For comparison we also show the correspond- 
ing result for gapped graphene, with a mass term equal 
to half the band gap of the GAL, A ~ 0.58 eV. In 
both cases, we fix the chemical potential at the lower 
band gap edge, i.e., fj. = A. We find reasonable agree- 
ment between the full GAL calculations and the simpler 
gapped graphene model. However, we note that a dis- 
tinct feature of the GAL structure is the additional reso- 



nance near Hui = 1.65 eV, which is absent in the simpler 
gapped graphene Dirac model. This resonance occurs 
due to transitions between bands that are not present in 
a simple two-band gapped graphene model of GALs. We 
will explore the details of the magneto-optical response 
of graphene antidot lattices in future work, and include 
the result here mainly to emphasize the power of the per- 
turbative formalism presented. 



V. SUMMARY 

A perturbative approach to calculating the optical Hall 
conductivity of graphene structures has been presented 
and applied to tight-binding models of graphene and 
graphene antidot lattices. The optical Hall response of 
graphene shows significant deviations from a simple Dirac 
treatment. While the Dirac model predicts a Hall con- 
ductivity of identically zero for a chemical potential at 
the Dirac point energy, results from our next-nearest 
neighbor tight-binding model indicate clear resonances at 
the band gap. The numerical results suggest that these 
effects should be measurable in experiments. Results for 
graphene antidot lattices illustrate that in this case, even 
the simple nearest-neighbor tight-binding model gives 
qualitatively different results than a simple Dirac approx- 
imation. 
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Appendix A: Linearized graphene 

Linearizing the tight-binding Hamiltonian of gapped 
graphene near the K point results in the celebrated Dirac 
approximation of graphene, 



H = 



A 3/2 (Kx — iK-y 

3/2 + iny) —A 



(Al) 



where we have introduced Ki = thiac, with Oc the nearest- 
neighbor distance = o/a/S We use this form of the 
Hamiltonian to evaluate the trace, noting that the lin- 



earization means that Tx 



(2) 



0. In polar coordinates, 



{Kx,Ky) = k(cos0, sin0) we find, after integrating over 



7 



the angular component, 



'Tr[T^y(z)] = 
10368TTt^ain{n - 2z) [z{n - z) + A2] 



9k2-4(z2_ A2) 



(A2) 



where Q = huj is the photon energy. We now use the 
residue theorem to perform the contour integral over z, 
yielding 



dz if{z) dc/) Tr [Txy{z) + T^y{z + O)] 



217V9k2 + 4A2 (9k2 + 4A2 - 172) 



/' (-VS^VITA^) - /' (v/9k2/4 + A2)](A3) 



where poles at z = 217 ± ^9^2/4 + A2 = 217 ± e(fc) have 
been ignored because, as discussed in the paper, the con- 
tour explicitly excludes these points. Inserting this result 
in Eq. (1) of the paper and converting the Brillouin zone 



integration to an integral over energy, this leads to Eq. (|4]) 
of the main text. Taking as their starting point the Lan- 
dau level structure of gapped graphene, Gusynin et al. 
have previously derived the off-diagonal magneto-optical 
conductivity of gapped graphenei^ Their result is stated 
as a sum over Landau levels, 



2n2^2 



J2 ([/(-^n+i) - fi-en)] - - f{en+l)]) 



n=0 



1 - 



A2 



A2 \ 1 



172 



(A4) 



where the energies are e„ = + nh?uj'^ for n > 0, 

yielding e„+i = e„^l + ^ e„ + h^u)l/{2en) in 

the low- field limit. Thus, in the continuum limit ^ = 
;i2c^2/2. Replacing /(e„) - /(e„+i) ^ -/'(cn)^ and 
converting the sum to an integral via ^ ~^ ! 
recover Eq. Q. 
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